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We explore the requirements for a Lyman-a forest survey designed to measure the angular diam- 
eter distance and Hubble parameter at 2 < z < 4 using the standard ruler provided by baryonic 
acoustic oscillations (BAO). The goal would be to obtain a high enough density of sources to probe 
the three-dimensional density field on the scale of the BAO feature. A percent-level measurement 
in this redshift range can almost double the Dark Energy Task Force Figure of Merit, relative to 
\& ' the case with only a similar precision measurement at z ~ 1, if the Universe is not assumed to be 

flat. This improvement is greater than the one obtained by doubling the size of the z ~ 1 survey, 
with Planck and a weak SDSS-like z — 0.3 BAO measurement assumed in each case. Galaxy BAO 
surveys at z ~ 1 may be able to make an effective Lya forest measurement simultaneously at min- 
imal added cost, because the required number density of quasars is relatively small. We discuss 
the constraining power as a function of area, magnitude limit (density of quasars), resolution, and 
signal-to-noise of the spectra. For example, a survey covering 2000 sq. deg. and achieving S/N = 1.8 
t | per A at g = 23 (~ 40 quasars per sq. deg.) with an R > 250 spectrograph is sufficient to measure 

both the radial and transverse oscillation scales to 1.4% from the Lya forest (or better, if fainter 
magnitudes and possibly Lyman-break galaxies can be used). At fixed integration time and in the 
sky-noise-dominated limit, a wider, noisier survey is generally more efficient; the only fundamental 
upper limit on noise being the need to identify a quasar and find a redshift. Because the Lya forest 
is much closer to linear and generally better understood than galaxies, systematic errors are even 
less likely to be a problem. 
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I I. INTRODUCTION 

o 

H 

The existence of new physics that causes the Universe to accelerate at late times is well established, both by 
observations of Type la supernovas d 0] an< ^ by combinations of other observables Q . The focus now is on probing 
the properties of the acceleration, most commonly thought to be caused by dark energy, a substance we know almost 
• i-h , nothing about except that it must have negative pressure. We often parameterize dark energy by the equation of 
■ state, w — p/p, where p is the pressure and p the density (w = — 1 for a cosmological constant). In general, w can be 
H | time-dependent 

Acoustic oscillations before recombination lead to a feature in the matter correlation function at the scale of the 
sound horizon at decoupling CMB observations pin this scale at 143 ±4 Mpc Recently, ^(| detected the 

expected enhancement (at > 3cr) in the correlation of the Sloan Digital Sky Survey (SDSS) Luminous Red Galaxies 
(LRGs). In coming years, precision measurements of the BAO feature will be used as a standard ruler to probe the 
equation of state of dark energy [IJ[ll[Il[ll[ll[li[llE 

We define the usual cosmological parameters in addition to w: The scale factor is a = 1/(1 + z), with Hubble 
parameter H(z) = a/a with H(Q) = Hq = h 100 km s _1 Mpc _1 . The baryon density is u>b — SI&/1 2 where fif, = pb/p c 
is the fraction of the critical density, p c = 3Hq /8ttG, in baryons. The matter density is uj m — fl m h 2 , where fl m is the 
fraction of the critical density in baryons plus cold dark matter. The fraction of the critical density in dark energy is 
£l w . The curvature is parameterized by 0^ = 1 — tt m — Q w . The primordial power spectrum is parameterized by its 
amplitude A and power law index, n s (we will not use tensors or running of the spectral index in this paper). The 
Thompson scattering optical depth to the CMB surface of last scattering is r. We will sometimes use 9 S at the CMB 
last scattering redshift, defined below in Eq. (J2J, as a free parameter in place of h. We measure 9 S in degrees so a 
typical value is 0.6. 
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In general, the quantities measured most directly in a BAO study are 

v s (a) = H{a) s a (1) 

in the radial direction and 

Daw 

in the transverse direction, where s is the comoving sound horizon at decoupling |25|. and Da(cl) is the angular 
diameter distance. For flat models, Da(cl)/cl = r{a) and 9 s (a) = s/r(a) with 

Ja a' 2 H(a') 

For curved models, 

r{a) 



(1 + z) D A (a) = R sinh 
with 



R 



(4) 



Finally, 



R = -= . (5) 



[** c s (a) 
J a z H{a) 



where a r ~ 1/1090 is the recombination redshift ( 25] gives a useful fitting formula for a r 's relatively small parameter 
dependence), and c s (a) — c/y/3 (1 + 3/9f,/4p 7 ), where p 7 is the photon density. 

In this paper, we discuss the possibility of measuring the BAO feature in a Lya forest survey (see also |2(|). The 
acoustic scale 143 Mpc corresponds to ~ 1.3 degrees on the sky at z ~ 2.8, our typical mean redshift. On this large 
scale, the clustering of the absorption can be assumed to follow a linear bias model in a way analogous to galaxies. 
This has been studied in simulations by |27| . The primary difference from galaxy surveys is in how the density field 
is sampled on relatively small scales. With galaxies, we observe a collection of discrete points, which we believe are 
a Poisson sampling of the underlying density field. For the Lya forest, we observe a collection of skewers through 
the density field, obtaining a detailed picture along the skewers, but no information in between. The reader who is 
uncomfortable with this should consider that they are probably perfectly comfortable constructing a density field from 
the ultimately zero dimensional galaxy sampling - in the forest we not only have a fully one-dimensional field, we even 
measure a continuous value rather than the discrete presence or absence of a galaxy. For galaxies there is intrinsic 
noise due to the Poisson sampling, while for the Lya forest there is intrinsic aliasing-like noise due to the discrete 
transverse sampling. As we will show, if the sampling of the Lya forest was only zero dimensional, this noise would 
take approximately the same white (uncorrelated) form as the galaxy noise. Of course, the Lya forest also has the 
advantage of being sensitive to the near- mean-density intergalactic medium (IGM) , rather than the highly non- linear 
positions of galaxies, so we can more confidently back up our assumptions with near-first-principles calculations. This 
is not to say that the Lya forest is decisively superior to galaxies as a high-z BAO probe. Ultimately, both should 
work, and the resource requirements to conduct an equivalent survey are the most important factor in deciding which 
is better. 

Computing the expected errors obtainable from the Lya forest as a function of survey parameters is the main aim 
of this paper. In SjH]we describe our basic method for estimating obtainable errors on 9 S and v s . Then in illlll we give 
the results. Finally, in fllVI we discuss the usefulness for constraining dark energy and curvature of a measurement of 
9 s (z ~ 3) and v s (z ~ 3). 



II. SETTING UP THE CALCULATION 



In this section we discuss our method for estimating the constraining power of future surveys. First, in W Al we 
define the parameters that describe the survey configuration. Then, in £111 Bl we explain our Fisher matrix calculations. 
Finally, in £111 CI we explain our theoretical model for the Lya forest power spectrum. 
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A. Assumed data set 



We assume a square survey with area A and comoving number density of quasars n q . We usually ignore evolution 
across the redshift extent of the survey, e.g., n q is the density at the central redshift. One would of course need to 
consider evolution in an analysis of real data, but, as we will show, the resulting error bars will not be increased 
relative to a Fisher matrix estimate that ignores evolution. n q will be a function of the central redshift and magnitude 
limit of the survey, based on the luminosity function, -y^-(z), taken from j2^. Beyond setting n q , the luminosity 
function is relevant for setting the distribution of noise levels in spectra. We assume that the survey probes the 
forest in the observed wavelength range A m ; n < A < A max , which we translate to z m ; n and z max using A = A Q (1 + z) 
(with A Q = 1216A). The central redshift is just (z m in + z max )/2. We assume the usable rest wavelength range in a 
single spectrum is A rcst , m in < A rcst < A res t,niax, generally using A rest , m in = 104lA and A rcst , m ax = 1185A (this is very 
conservative - one could probably extend the baryonic oscillation analysis into the Lyman- /3-influenced region without 
much difficulty). We will occasionally call the length in the IGM probed by a single line of sight L q (~ 330 h~ x Mpc 
for these limits). The spectrograph is assumed to have rms resolution er^ (FWHM = 2.355 cjr = X/R). Spectral 
pixels are assumed to have full width l p = <jr. We assume rms noise per pixel ct/v, in units of the mean transmitted 
flux level in the forest. 



B. Fisher matrix calculation 



We estimate the errors on parameters pi, obtainable using a future survey, to be (F 1 )u, where Fij is the Fisher 
matrix [2^| : 

with 

2£ = lndet C + S T C~ 1 S + constant (8) 

for a real mean zero Gaussian vector S of data points with covariance matrix Cy(p) = (SiSj). From [2jj (and references 
therein), we have 



F i3 ^-TrfC" 1 — C- 1 — ) . (9) 



We could always take Si to be individual pixels in spectra and evaluate Eq. by brute force; however, this calculation 
would quickly become difficult. For example, our fiducial survey would have about 40 million pixels. While we could 
probably reduce this number enough to produce a reasonably invertible covariance matrix by scaling up from a much 
smaller survey and pushing the pixel size and angular sampling to just barely resolve the BAO feature, it is better to 
look for a more efficient method. 

Methods for applying the Fisher matrix formalism to a BAO survey using galaxies are well-developed . Counts 
of galaxies (say, in cells smaller than the scale of interest) are the density field of interest, but the calculations are 
done in Fourier space. Assuming an approximately uniform selection function, the covariance matrix of Fourier 
modes is approximately diagonal if the wavenumber kj is discretized in bands of width 2?r/ Lj, where Lj is the 
survey width in direction j. The covariance matrix of the real and imaginary parts of Fourier modes is simply 
Cii(p) = \ [Pg0ti,p) + Pn], where P g is the redshift-space galaxy power spectrum and Pjv = 1/h s is the Poisson 
noise power for mean galaxy density n g . Simple formulas for the Fisher matrix as an integral over k can be derived. 

The Lya forest is slightly more complicated. Here, the quasar lines of sight going through some volume are 
approximately random. Our density estimate for the volume is not determined by the number of lines of sight probing 
it, but rather by some average of the absorption in those lines of sight. Working toward a covariance matrix of Fourier 
modes, we will take Si in the Fisher matrix to be the Fourier modes of the weighted density field 

*(x) = ^[i f (x)+i ff (x)]. (10) 
w 

where Sp{x) = Sp(z,8) — F(z,9)/F — 1 where F = exp(— r) is the transmitted flux fraction in the forest, Sjy(x.) 
is the spectral noise, u>(x) is the weight at x, and w is the mean weight. We will not specify the weights at this 
point, but generally they should depend on pixel noise variance. The relevant quantities are all defined over fully 
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three-dimensional space - if some location is not probed by a quasar line of sight, the weight is simply zero. The 
correlation of Fourier modes of S(x) is 

(Wk'} = /cbcdx'exp(£k-x + £k'-x') [£ F (x-x')+6Kx,x')][l + Mx)][l+Mx')] » ( n ) 

defining 5 w (x) = w(x)/w - 1, and (, x (x,x') = (Sx(x)5 x (x')) . 

We do not know the weights in advance, because they depend on the random quasar locations and luminosities. 
We can, however, evaluate the average of the weight term, which should give us a good estimate of the constraining 
power of the survey (for a typical survey of ~ 10000 quasars, a quantity measured from a single realization should 
not deviate much from the average quantity). Our pixel weights will be approximately uncorrelated in the transverse 
directions, and perfectly correlated in the radial direction (for simplicity we are assuming pixels in the same spectrum 
have the same noise). For a continuous white noise field, £(x) = P S D (x), where P is the constant power spectrum 
of the field. This gives 

(Mx)Mx')) =p™s d (x±) (12) 

where P 2B is the power spectrum of 8 W in the two angular directions. This means 

(6 k 5 k ,) ~ (27r) 3 ^(k + k') [P F (k) + ^ D (fc||)P^ D + P C N S ] , (13) 

where P^ D (fc||) is the usual ID flux power spectrum measured along single lines of sight, and P^f is the weighted 
noise power, which does not separate like the other term because the weights depend on the noise amplitude. Note 
that modes with different k are still approximately uncorrelated. 

We now need to evaluate the P 2B term. The low-k power spectrum for a pixelated white noise field is P = <r 2 v, where 
a 2 is the pixel variance and v is the pixel volume (in a general sense, i.e., here the volume is actually an area because 
we are working in two dimensions). As a computational device, let us imagine pixelating the transverse directions 
into small cells of angular width I (we will eventually take I — ► 0). Then, P 2B = cr^l 2 , with a 2 ^ = (u> 2 ) jw 2 — 1. The 
noise in a spectrum is a function of the apparent magnitude of the quasar, so the mean weight is 



drJ^- L q I 2 w(m) = l 2 L q h . (14) 
dm 



Note that L q l 2 dn q / dm is the probability that a given small volume of the IGM will be probed by a line of sight to a 
quasar of apparent magnitude m. Similarly, 



Therefore, in the limit I — > 0, 



A similar calculation leads to 



where 



(w 2 ) = I dm^pl L q I 2 w 2 (m) = l 2 L q I 2 . (15) 
' , dm 



p™ - (is) 



Pf = 4^ (17) 



dm^- a 2 N (m)w 2 (m) . (18) 
, dm 

Note that l\L q jl v is something like an effective 3D pixel density (recall that l p is the pixel width), while I3/I1 is the 
effective noise variance in these pixels. 

We take the weights to have the simple Feldman, Kaiser, & Peacock (hereafter FKP) inverse variance form |30l |. 



(19) 
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where P$ is the typical signal power and P/v(x) = Pjv[rn(x)] is the noise power level associated with the noise level 
at x (if there is no quasar probing x, this level is infinite so the weight is zero). The numerator is chosen to make 
w{x) = 1 in the low noise limit. There is a subtlety at this point in that P/v(m) is not unambiguously defined. If 
quasars were evenly distributed in angle, and all had the same magnitude, then clearly P/v = aJfl p /n q L q . The only 
reasonable possibility in the realistic case seems to be to take Pjv(m) = a^{m)l p / I\L q = crj f (m)/n p s , where is 
the effective pixel density. Note the similarity to galaxy shot noise: if the overall mean density of galaxies is n g , and 
the mean at point x is n s (x), then the Poisson- noise variance at point x is Up — n g /n s (x) for cells with volume 
1/rig, and the noise factor in the usual FKP weights is Op{x)/n g (recall that our m dependence is equivalent to x 
dependence). Since n° ff and the weights are mutually dependent, we determine them using a few iterations. After 
some experimentation showing that the results are not very sensitive to the choice of P<j, we take the shortcut of 
using one constant for P$: the total flux power in the central model at k = 0.07 h Mpc^ 1 , fi = k\\/k = 0.5, including 
the aliasing term (which is formally signal, although in practice it acts as noise). 

There are a couple of loose ends in this derivation. First, we are going to assume that, having computed the 
power spectrum, i.e., the covariance matrix of Fourier modes, averaged over realizations of the weights, we can simply 
proceed with the calculation of the constraining power of a survey as if this were the exact power spectrum; however, 
the averages of the covariance matrix of power measurements, or of the Fisher matrix, are not necessarily the same 
as what one computes from the averaged power spectrum. Ideally we would only average at the last possible point in 
the computation (remember, the Fisher matrix is already an average of the derivatives of the likelihood function over 
all possible future data sets). Unfortunately, without the assumption that the Fourier modes are uncorrelated, we 
are stuck with the problem of inverting a very large covariance matrix, and now we would need to do it many times 
to average over realizations. To investigate this problem, we computed the average covariance matrix of the power 
spectrum measurements and found that the new terms are suppressed by a factor of the number of quasars in the 
survey (i.e., the ratio of the area of the survey to the area per line-of-sight) relative to the standard terms, as long as 
the second term in Eq. I|13(l is not much larger than the first. This will be the case any time the BAO measurement 
is possible, so it seems we are safe to ignore this issue. Second, with no Nyquist frequency for a randomly sampled 
survey, it is not clear at what k we should stop counting modes. This is worrisome because we will find that our 
constraining power will continue to increase somewhat beyond the transverse Nyquist frequency associated with the 
effective mean density of lines of sight. This means that we are using more Fourier modes than we have data points, 
and the assumption that they are uncorrelated probably has to break down somehow, e.g., due to an accumulation 
of small correlations. To be conservative, we only use modes with k less than the Nyquist frequency associated with 
the effective mean density. We also eliminate modes with k > 0.5 /iMpc -1 (relevant only in the radial direction). To 
allow for continuum fitting, we drop all the modes having the first 27Vq.i os discrete values of k\\ , counting both positive 
and negative k\\ , where i\r qi i os is the number of spectra it takes to cover the redshift range of our survey (the first 5 
values of ku for our standard survey). These last two cuts make no noticeable difference in the results. 

We note for the future that one would probably perform the data analysis using a more sophisticated maximum 
likelihood method [29j, but this FKP-like method should be good enough for the Fisher matrix error estimate. 

C. Three-dimensional Lya forest power spectrum 

Missing from the discussion so far has been any mention of how we predict the Lya forest power spectrum and 
its parameter dependence that is needed for the Fisher matrix calculation. Fortunately, this is basically a solved 
problem: We use the fitting formula for Pp(k, p) from j^j. We start with a standard ACDM transfer function from 
CMBfast H3 with a$ = 0.897, fl m = 0.281, tt b = 0.0462, h = 0.710, and n = 0.980 The amplitude and slope of 
the primordial power spectrum (i.e., a$ and n) are free parameters to be marginalized over, while we leave the other 
parameters fixed because, from the point of view of the Lya forest, they are either well constrained (f2j,) or degenerate 
with the slope and amplitude. Additionally, we marginalize over four parameters of the Lya forest model: the mean 
transmitted flux level, F, the temperature at the mean density, T1.4, the slope of the power law temperature-density 
relation, 7 — 1, and the large scale anisotropy parameter (3 |22|. These parameters set the large scale bias of the Lya 
forest power spectrum (the bias is effectively free because it is hyper-sensitive to F). (3 is given by |53 as a function 
of the other parameters, but we choose to marginalize over it to be sure we are not relying on this prediction. Finally, 
the parameters of interest are the radial and angular scale factors, v s (a) and 6 s (a), which we use to convert the linear 
power spectrum from comoving Mpc/h units to the observed velocity (redshift) and angular coordinates. There is 
little chance that our power spectrum predictions can b e sig nificantly wrong, because they agree with high precision 
measurements of the ID power in single quasar spectra |33l l34l l35l I36L l37j . Note in particular that the aliasing- like 
noise term in Eq. (|13|) . which turns out to be critical to our calculation, is precisely the well-measured ID power. 
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III. RESULTS 

There are several interacting degrees of freedom for a survey so to give us a concrete start we imagine the Lya forest 
survey piggy-backed onto something like the proposed WFMOS survey of low redshift galaxies [2(|. This proposal 
is to observe two million galaxies in the redshift range 0.5 < z < 1.3 in a 2000 sq. deg. area using an 8m telescope 
with an R=2000 spectrograph with exposures of 30 minutes. We take the number density of quasars as a function of 
luminosity from |2g, finding that g magnitude limits (21, 22, 23, 24, 25) correspond to (8, 20, 41, 77, 136) quasars per 
sq. deg. within the relevant redshift range. We assume the noise is sky-dominated and for our baseline WFMOS-like 
survey assume S/N=(ll, 4.5, 1.8, 0.7, 0.3) per lA for g=(21, 22, 23, 24, 25). There are of course a substantial number 
of quasars brighter than any given magnitude limit. We track the S/N over the source counts, rather than assuming 
that all quasars are at the limiting magnitude. At g > 23, the number density of Lyman-break galaxies (LBGs) 
starts to become high enough to provide significant extra sources. In cases where we include these, we take the LBG 
luminosity function from [38|], assuming g=R+0.7, finding (0.3, 116, 2325) LBGs per sq. deg. for g < (23,24,25) 
(m^ = 24.54 for these galaxies). 

Figure n shows the first results, where we assume that we are observing the Lya forest in the wavelength range 
3900-5229A (2.2 < z < 3.3). The WFMOS galaxy survey is expected to obtain la errors of 1.0 and 1.2% on 
Da(z) and H{z), respectively, at z ~ 1, and, more relevantly, 1.5 and 1.8% at z ~ 3 from a separate, similarly costly, 
high z galaxy survey reaching R < 24.5. Our errors depend on the magnitude limit, reaching better than 1.4% in both 
directions for g < 23. The results improve more slowly with increasing magnitude limit because the faint spectra, at 
fixed observing time, are becoming too noisy to be very useful. If very faint quasars and LBGs can be identified, the 
errors could be as small as 0.5%. Note that the Lya forest BAO measurement generally constrains H(z) better than 
Da{z), the opposite trend from galaxies |20|. 

To make the results clearer, we show the obtainable constraints on band power measurements in Fig. [21 (our standard 
constraints on the BAO scale do not go through this intermediate step). We see how the inclusion of aliasing- like 
power from the discrete transverse sampling, i.e., the 2nd term in Eq. I|13l) . dilutes the amplitude of the BAO 
features, especially for modes transverse to the line of sight which are not enhanced by large-scale peculiar velocities. 
The absence of low-k transverse modes can be traced, somewhat counter-intuitively, to the fairly large minimum fc|| 
resulting from the relatively small radial extent of the survey compounded by the loss of modes dropped to allow for 
continuum fitting. The absence of high k transverse modes is due to the sparse transverse sampling. Spectral noise 
power is, by construction, always subdominant (~ 50% of the signal power including aliasing power), because quasars 
faint enough to contribute a lot of noise power are discarded by the weighting. 

For comparison, in Fig. [3] we imagine piggy-backing the Lya forest survey onto the high-z galaxy survey of |20| . 
which is proposed to cover 300 sq. deg. with 240 minute exposures, obtaining errors of 1.8% on H{z) and 1.5% on 
Da (z) using galaxies (note that the S/N we assume is guided by but should not be taken as a prediction for WFMOS 
sensitivity). We find similar errors to those from galaxies at a similar limiting magnitude, but the wider, shallower 
survey is probably a better option for the Lya forest. These figures uncover an important issue for survey planning: 
for a fixed survey configuration, the results are quite sensitive to the limiting magnitude, i.e., how bright does a quasar 
have to be to be identifiable? 

We now explore the optimization of the survey. The errors scale precisely as A -1 / 2 down to surveys much smaller 
than the ones we are discussing. It seems likely that any spectrograph used for this project will have more fibers than 
necessary, so the main question is how long to integrate before moving on to gather more area. Figure 0] shows the 
expected errors as a function of the S/N obtained per A for a g=22.5 quasar, with a corresponding rescaling of the 
survey area to keep A (S/N) 2 fixed. As in previous figures, we are including a realistic distribution of magnitudes, 
i.e., the noise is higher for fainter quasars and lower for brighter ones. We see that a wider, noisier survey is generally 
superior. Limitations on this progression will be set by telescope overhead and the minimum S/N needed to identify 
a quasar. We isolate the change with noise level in Fig. \5\ 

In Fig. El we show the effect of decreasing resolution. Only relatively poor resolution is needed. R=250 is essentially 
as good as 2000, and 125 is not too bad. These numbers are easy to understand. At R=125, the power suppression 
factor at our FKP weight point, k = 0.07 /iMpc -1 , [i — 0.5, is only 11%. This suppression factor increases very 
quickly with fc|| ; however, the aliasing noise that we show in Fig. [21 is itself suppressed by limited resolution, so the 
effect is somewhat less than one might otherwise expect. 

Finally, we perform a few tests to make sure our calculation is robust. First, to test our code, we recompute the 
errors for the galaxy survey in j20j. We can do this easily by dropping the aliasing term, setting the noise variance 
in cells defined by the pixel size and mean separation of the spectra to the inverse of the mean number of galaxies 
in a cell, and replacing the Lya forest flux power spectrum with the galaxy power spectrum. We find errors 1.5% on 
Da{z) and 1.6% on H(z), assuming they used fi m = 0.35, as = 0.9, b = 3.3, h — 0.65, n = 1.0, following [ljfl. This 
compares well with [2(|'s 1-5% and 1.8%, respectively. 

One approximation we have made is to ignore evolution, i.e., to use the power spectrum, quasar number density, 
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FIG. 1: Upper thick lines show constraints as a function of g magnitude limit on the radial and transverse BAO scales for a 
survey similar to the proposed WFMOS low-z galaxy survey [201 ] . assuming 2000 sq. deg. and R=2000. Magnitude limits (21, 
22, 23, 24, 25) correspond to (8, 20, 41, 77, 136) quasars per sq. deg. for the |28| luminosity function, and we assume S/N=(ll, 
4.5, 1.8, 0.7, 0.3) per A. Lower thick curves add spectra from LBGs. Thin lines show the completely unrealistic case where we 
ignore the aliasing-like noise power caused by discrete sampling. 



etc. from the center of the redshift interval for the full interval. As a test, we split the sample into two redshift bins, 
and compute the Fisher matrix separately for each, including marginalizing over two independent sets of nuisance 
parameters. The errors on each bin increase of course, but the combination of the two bins actually gives slightly 
smaller errors than the original full Fisher matrix, because of small differences in the averaging. Not surprisingly, the 
low redshift half of the survey gives significantly smaller error bars than the high redshift half (the quasar density is 
higher at fixed apparent magnitude at low z). 
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0.1 0.2 0.3 

k [h/Mpc] 

FIG. 2: Error bars show the fractional error on the power in bands of Afc = 0.02 /tMpc" 1 , for the g < 25 case from Fig. Q 
(without LRGs). Lines show the ratio of Lya forest flux power for Qb ~ 0.0462 to Qt = 0.001. Black (solid line, error bars 
shifted slightly right) shows /i = k\\/k > 2/3, red (dashed) shows 1/3 < fj. < 2/3, and green (dotted, error bars shifted left) 
shows < 1/3. For comparison, the thin line shows the ratio of power without aliasing noise (renormalized for clarity). 

To test that we are really measuring the BAO feature and not some broadband feature in the power spectrum, we 
run our error computation using a transfer function with fl b = 0.001. We find that the errors for the survey in Fig. G] 
are always greater than 5%, i.e., our measurement is clearly based on the baryonic feature. Finally, we note that our 
results are completely insensitive to removing the marginalization over (3 or adding a marginalization over the noise 
amplitude. 

A 30 sq. deg. pilot study reaching g < 21 with S/N = 11 at g = 21 should detect baryon oscillations at ~ 2a, in the 
sense that the Fisher matrix prediction for a measurement of Hi,, marginalized over the nuisance parameters discussed 
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FIG. 3: Similar to Fig. Q except covering 300 sq. deg. with longer exposures, following the high-z galaxy survey proposed in 
ppf . g magnitude limits (21, 22, 23, 24, 25) correspond to S/N=(32, 13, 5.1, 2.0, 0.8) per A. 

above, predicts an error 0.023 for = 0.046. Reaching g < 22 with S/N > 4.5 would produce a 2.7a detection. A 
detection of this kind would not be fundamentally interesting; rather, we give these numbers to indicate the type of 
survey that would be needed to produce a good solid measurement of the Lya forest power on the appropriate scale. 
As usual, a survey four times larger with half this S/N would produce a better measurement. 
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FIG. 4: Similar to Fig.Q but investigating the effect of the trade-off between noise and area. From top to bottom (black, blue, 
cyan, green, magenta, red) we show, at g—22.5, S/N=(5.7, 4.1, 2.9, 2.0, 1.4, 1.0) per A. The results improve with decreasing 
S/N because in each case the area of the survey is A = 2000 (2.9 N/S) 2 sq. deg., assuming sky-dominated noise. Errors scale 
as A' 1 / 2 . 



IV. CONSTRAINTS ON DARK ENERGY AND CURVATURE 



In 3IV Al we give projected constraints on the most standard specific parameterization of the dark energy. In 3IVBI 
we discuss the general usefulness of a BAO measurement at z > 2, in a non-parametric way. 
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: — d a (z) s -> ; 

H(z) s 

I i i i i I i i i i I i i i i I i i i ■ 

21 22 23 24 25 

§max 

FIG. 5: Similar to Fig. Q but investigating the effect of the changing the noise level at fixed 2000 sq. deg. survey area. From 
bottom to top (black, blue, cyan, green, magenta, red, black) we show, at g=22.5, S/N=(5.7, 4.1, 2.9, 2.0, 1.4, 1.0, 0.72) per A. 

A. Parametric models 

We discuss first the usefulness of a BAO measurement at z ~ 2.8 for constraining dark energy with equation of 
state w(a) — wq.6 + (0.6 — a) w\ with wqs and w' as parameters. We saw above that the true center of weight of the 
Lya forest survey would probably be a bit lower than z ~ 2.8 but we find negligible sensitivity of the hnal constraints 
to the exact Lya forest redshift. The pivot point a = 0.6 was chosen to make w a and w' roughly uncorrelated for 
our scenarios. We use the figure of merit (FoM) defined by the Dark Energy Task Force ^3 (DETF) as our primary 
measure of survey value. The DETF FoM is simply the inverse of the area within the 95% confidence contour in the 
w a — w' plane, which has the important characteristic of being independent of the chosen pivot point. 



12 




FIG. 6: Similar to Fig. Q but investigating the effect of changing resolution. From bottom to top (black, blue, green, red) 
lines show R=(2000, 250, 125, 62.5). 



We include the projected constraints from the Planck CMB experiment. We use standard Fisher matrix techni ques 
[3^. |40| , following ^jj in ignoring foregrounds but using only the 143 GHz channel. Our results agree well with |4l| 
when we use similar parameter combinations. For reproducibility, our results are given in Table [I] We see that the 
CMB alone only weakly constrains dark energy (note here that the Fisher matrix technique may not give reliable 
errors in cases where the errors are large, but this should not be a problem once other data constrains these degenerate 
directions). 

Table ITT1 shows the FoM and constraints on wo.e, w' , and sometimes fife, for several different combinations of data. 
To avoid sensitivity to broad-band power, we use separate (independent) slope and amplitude parameters for the 
CMB and Lya forest. We use three different BAO constraints: First, we always assume that SDSS will produce a 
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TABLE I: Assumed values, errors, and error correlations from Planck. The first row is the assumed value of the parameter, 
the second is the error, and the rest are the correlation matrix. Very large errors should be interpreted only qualitatively. 
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0.034 


0.577 
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0.161 


-0.117 
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0.180 
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0.315 
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-0.645 
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-0.004 
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0.069 
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-0.010 


0.684 
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0.252 
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-0.059 
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-0.018 
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-0.006 


0.315 


-0.019 
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TABLE II: Errors on wo.6, w' , and fit [with w(a) — rao.6 + (0.6 — a) w'] for different data combinations, along with the DETF 
figure of merit (FoM), defined to be the inverse area inside the 2 — a contours in the too. 6 — w' plane. BAO errors are 5.8% 
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text). 
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0.0038 
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0.52 
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1.40 
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0.27 
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Y 
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1.31 


0.074 


0.87 
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Y 
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Y 


0.54 


0.122 


1.35 
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Y 


Y 


Y 


1.61 


0.067 


0.74 





5.8% error on the radial scale at z — 0.3, v s (z — 0.3), and 5.2% error on the transverse scale, 9 s (z = 0.3) |17j . We 
optionally add a 1% constraint on both distance scales at z = 2.8. Finally, we optionally include a ~ 1% constraint 
at z ~ 1. To avoid artificially degrading the value of a z ~ 1 survey, we spread this measurement over points at 
z = 0.6,0.8,1.0, and 1.2, using the relative errors from p"?| but modifying the overall normalization to make the 
combined error from the eight points 0.7%, i.e., the same total precision as the two 1% z = 2.8 measurements. We 
choose a simple 1% because different survey configurations can lead to many different combinations of errors. The 
point here is primarily to study the general usefulness of constraints at different redshifts, not any specific survey. 

We see, as expected, that in the presence of Planck constraints the measurement at z = 2.8 is less valuable for 
constraining dark energy than the measurement at z ~ 1; however, if £lk is allowed to vary, adding the z = 2.8 
measurement nearly doubles the DETF FoM. This improvement is larger than the affect of doubling the size of the 
z ~ 1 survey. The higher z measurement is actually slightly more valuable than z ~ 1 for constraining Q}-. If we 
assume flatness the improvement is more modest, although not completely negligible. If the CMB is taken out of play 
for some reason, the z = 2.8 measurement becomes more valuable than z — 1, although in that case both are really 
needed to provide an interesting measurement. 

We note that the commonly used approximation that the Planck constraint can be represented by the constraints 
on 9 S (0.04%, or perfectly known) and w m (1%) alone, treated as independent of each other and other parameters, 
works very well, in the sense that the resulting FoM agrees with the full-Fisher matrix version to better than 10%. 

B. Non-parametric 

We will next consider the measurement of the acoustic scale at z > 2 in the context of testing the fiat ACDM 
model. We will expand in small perturbations around the fiducial model, working to lowest order in the non-constant 
dark energy and curvature. Hence, we write H 2 = (87rG/3)(pfid + px) + QkHq (1 + z) 2 , where pad is the sum of the 
matter, radiation, and cosmological constant contributions and px includes the dark energy density that differs from 
the z = value. We want to test whether px = or not, keeping in mind the uncertainties in pad and curvature. 
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1. Transverse Scale 



In the transverse direction, we measure i z = 9 s {z)^ 1 = (1 + z)Da(z)/ s. For example, we measure ^ios9 to wonderful 
accuracy (0.35% with 3-year WMAP, <C 0.1% with Planck). Now we consider a measurement at z w 2 and construct 
£ios9 — £ z as a way to isolate the new information beyond that available in the CMB. To lowest order in curvature, 
the angular diameter distance 



(1 + z)D A = r + 



(20) 



This then yields 

£i = 1090^(1089) - (1 + z)D A (z) ^ K1089) - r(z) + Ogg _ ^ 

s s 6sc z 
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1089 cdz Qgflg 
-ff 6sc 2 
Q K £T 2 r 2 (1089) 
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(1 + d + rf 2 ) 
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ff fid 2s 



10S9 



cdz 

cdz ( p x 3fl K H 2 (l + zf 



Hfid \Pfid 



8nGp 



'fid 



(21) 



where we've define d — r(z)/r(1089) and assumed that \CIk\ "C 1 and <C pad (suitable for z > 2). The latter 
means that px/pad = fix- Keeping only lowest order in the perturbations around the fiducial model, we also have 



3fl K H 2 (l + z) 



fi K 



8nGp 



lid 



n m {i + z) 3 + n A n m (i + z ) 



[l-2A + 0(A 2 )] 



(22) 



where A = Cl\/2D. m (l + z) 3 . We can drop the higher orders in A when working at higher redshift. 
We now define a weighted average of fix 



fl x = 



1089 cdz 

s HflA 
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z Hfid 



(23) 



and rearrange to find 

fix = 2 



2s 



1089 cdz 
Hfid 



(4 



10S9 



K 



fr 2 (1089)(l 



d + d 2 ) - 



6A 
7 



m m {i + z) 



(24) 



The last term comes from doing the integral for the curvature term in (|21|l . For standard cosmologies, the coefficient 
in square brackets is about 5 for z — 2-3. The middle term must be —2 in the fiducial model, because the result 
cancels to fix = 0. 

The question is now how accurately we can measure fix in light of the uncertainties on the various terms on the 
right-hand side of the equation. In particular, we want to know how the uncertainties in fl m h 2 , fl m , and fl\ enter 
(i.e., how well specified the baseline model is). It is useful to rearrange the middle term as 



flm,H() 



1089 



c dz 



Had 



(£z — ^1089) 



(25) 



The combination s(ri m 7J 2 ) 1 / 4 is picked to cancel out most of the dependence of s on fl rn HQ. We can manipulate the 
integral as 
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Hfid ^ffl^H 2 
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dz 



^{i + zf + {fi r /fi m ){i + z y 



(1 + *)- 
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(26) 



If we expand the square bracket terms and neglect the radiation term when integrating the terms with A, then we 
find 



1089 



c dz 



Hfid yJfl m Hl 



1089 dz 2_ /A _ 3A^ 

^{l + zY + {fl r /fl m ){l + zY y/T+l\7 26 



(27) 
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2c 



where a = 1/(1 + z). With this, we find 



e 1 1 1089 



A 3A 2 



(28) 



tlx = 2 
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O m (l + z) 



(29) 



Now we can look at the errors in these terms. We denote cr(x) as the standard deviation of x. The error in £l m h 2 is 
expected to be below 2% with Planck data. That means that the fractional error in {n m H%) 1 / A is below 0.5%. The 
quantity s(ri m 7J 2 ) 1 / 4 essentially depends only on flbh 2 , but this is to the 1/8 power, so the error will be below 0.1%. 
For A, 



a(X) 



1 



a(Q r< 



2(l + z) 3 fi r 



(30) 



As we care about A/7, the value ends up being about a(X/7) — <r(f2 m )(l + z)~ 3 . With today's cosmological constraints 
( e -g-i @jEJ)j this error is about 0.1% at z — 2; of course, this will shrink in the future. In other words, our low 
rcdshift data constrains A well enough that the uncertainties in the extrapolation to z > 2 are tiny. The contribution 
of A in the curvature term is smaller yet; we will drop these. The errors in the radiation term z eq are very small. 

This leaves the fractional error in t z — 4o89- The error in will be very small with Planck, so we neglect it. 
Then 



a(£ z -hos9) _ <r{lz) 



1 



J.0.V1 



4 ^1089/4 -1 



(31) 



The quantity £wag/£ z — 1 is 1.7 for z = 2 and 1.2 for z = 3. 

Hence, the error in Clx is dominated by the quadrature sum of half the fractional error in £l m h 2 (< 1% for Planck), 
the fractional error in l z times 2/ (4o89/4 — 1) (which is 1-1.7, depending on redshift), and five times the uncertainty 
in the curvature Q.k- Of course, one might opt to move the curvature to the other side of the ledger and constrain 
Qx — 5Hx- 

It is a good approximation to think of Vtx as 



x 



f a din a a 1 



/2 



(32) 



One might compare this to the effect of anomalous dark energy on the growth of structure, which enters as a suppression 
of growth approximately as a factor 



I — — / d In a Clx 



(33) 



The weightings differ by a 



1/2 



2. Radial Scale 

Next, we turn to the radial acoustic scale. Here we are measuring v s (z) = sH(z)/(l + z) at a given z. Performing 
the same expansion, we have 



v s (z) 
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If we expand in powers of A, we get 

v s (z) = [s(fl m H 2 )^ 4 ] (1 + *) 1/2 (0» Hi) 1 ' i 
Solving for fix, we have 



1 + A- y + 0(A 3 ) 
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fi x fi 

~2~ + 2fl m (l + z) 
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(36) 

Again, the middle term must be 2 in the fiducial model. 

Looking at the error budget, again we expect to know {flmHl) 1 / 4 to better than 0.5%. The error in A is about 
a(fl m )/4 at z = 2, which is already about 0.7% today and should drop toward 0.1%. The coefficient of flu will be 
about unity for a survey at z f» 3. So the error in fix is about the quadrature sum of half the fractional error in 
il m Ho (somewhat below 1%), twice the fractional error in the acoustic scale v s (z), and the error on fix- 



3. Measurement Goals 



It is not clear what quantitative goal one wants to set for the measurement of fix or fix- One goal is to simply 
detect the cosmological constant at high redshift, i.e., to exclude fl x (z) = — fl\(z). When one averages this model 
over redshift, one gets fix = $l\/7VL m (l + z) 3 . This is only about 1% at z = 2. Hence, it is very hard to detect the 
absence of high-redshift dark energy using the transverse distance scale. However, it is fairly easy with the radial 
scale: fix = — fiA/fi m (l + z) 3 is about 10% at z = 2 and 4% at z — 3. A 1% measurement of H{z) at z = 2 should 
detect the dark energy at 3 to 5 <r, depending on whether the error bars on flmh 2 are 2% or somewhat better. 

Of course, extra dark energy can be detected too. The 1% measurement of H(z) would bound (at 3-5 a) the energy 
density of the dark energy at z = 2 to be within a factor of two of its low redshift value. 

The above statements assumes a flat cosmology, but the sensitivity of the transverse scale to curvature is 5 times 
larger than that of the radial scale. Hence, one can use the transverse acoustic scale to control the curvature to levels 
well smaller than the errors on the fix measurement from the radial acoustic scale; this assumes that fix < fix- 
Putting this another way, one can mix the transverse and radial scale to produce a curvature- independent measurement 
of a combination fix + fix/5, the coefficient being appropriate to the redshifts considered here, Many dark energy 
models will not yield zero for this combination; this would require that fix is larger at redshifts above the survey 
redshift and of the opposite sign. Hence, one can perform a reasonably generic search for deviations from the 
cosmological constant at the survey redshift. 

Whether this is interesting depends considerably on one's model for dark energy. The w = wo + w a ({ — a) model has 
the unfortunate property of demanding that there are no changes in dark energy at high redshift that aren't heralded 
with even larger changes at low redshifts. If one relaxes that assumption it is certainly the case that we don't know 
the dark energy density at z = 2 to 30% rms. 

Another question is whether there is a steady value of fix at high redshift, as tracker models would have |4^. Here, 
fix is reasonably close to fix, and one could reach a error bar of about 1%, assuming a flat cosmology. One may 
compare this to constraint available from measuring the growth of structure relative to the CMB. Here, one is limited 
by the optical depth to z = 1000. With Planck, this may be measured as well as 0.5%, so if one could measure the 
growth function to that level, one could measure a value of fix in the 0.1-0.2% range (l-cr). This method, however, 
is at least partially degenerate with the suppression caused by non-zero neutrino mass |43|]. A suppression of 0.5% 
corresponds to about 0.015 eV, already well into the region indicated by the atmospheric neutrino results. Hence, the 
degeneracy between neutrino mass will be an issue as one tries to push constraints on fix below 1%. 

Looking beyond the scope of the surveys discussed in this paper, the full-sky cosmic variance limit on the acoustic 
scale at these redshifts is superb. At this point, the uncertainty on fl m h 2 from CMB experiments such as Planck 
will produce sufficient uncertainty in the acoustic scale as to dominate the error budget on fix or fix- From a 
statistical point of view, such BAO surveys themselves could measure fi m /i 2 better than the CMB, thereby restoring 
the precision in the sound horizon, but the systematics in broad-band power measurement from low-redshift surveys 
will surely be worse than those in the CMB. Another option is to form a combination fix + fix that is independent 
of fl m h 2 . This combination, however, is not independent of curvature; essentially one is using the angular acoustic 
scale in the CMB to calibrate fl m h 2 . Nevertheless, this is a simple example of the idea that the BAO measurement 
of Da(z) and H(z) does produce an internal cross-check, namely that Da is an integral of H(z), that can be used to 
eliminate certain nuisance parameters. 
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V. CONCLUSIONS 

As a probe of dark energy, baryonic acoustic oscillation features have a significant advantage over standard candles 
such as supernovas in that we have a very clear theory describing them. Even the imperfection associated with our 
lack of understanding of exactly how galaxies populate dark matter halos does not present a significant problem, 
because the scale of halos is well separated from the scale of the BAO. The Lya forest extends this advantage in that 
we can perform something relatively close to a computation from first-principles. 

In this paper we have computed the statistical power that can be expected from a large, three-dimensional Lya 
forest survey probing the BAO scale. If the Universe is not assumed to be flat, a measurement of the BAO scale at 
z ~ 3 provides a interestingly large improvement in the constraints on dark energy, as quantified in 9IVI using the 
Dark Energy Task Force figure of merit. A range of combinations of area and magnitude limit could lead to a viable 
survey, as described in 91111 Resolution and signal-to-noise ratio requirements are modest. It will probably be most 
efficient to piggy-back Lya forest surveys on top of lower redshift galaxy surveys. 

As with other probes, the Lya forest version of the BAO measurement should be robust against systematic errors. 
Long-range effects in the UV ionizing background will not produce preferred scales and can't mimic the acoustic 
peak. The poorly modeled highest density absorption (e.g., DLAs) will not produce a BAO signal cither. Continuum 
features could possibly create a bump on the appropriate scale, but these should be easy to control because they 
would be associated with fixed quasar rest wavelength ranges, unlike the BAO feature. In the worst case we would 
lose only a small fraction of pixel-pairs by completely ignoring correlation between pixels in the same spectrum. Since 
we would be using multiobject spectrograph techniques, fluxing errors would create density variations as a function 
of redshift for all objects in a region. This creates fake large-scale power; however, one can project this purely radial 
power out with minimal information loss |44l |45| . 

The ultimate Lya forest survey could include faint Lyman break galaxies to obtain a higher density of probes. 
However, getting a redshift for these galaxies is not trivial. For the 75% that don't have Lya emission, one needs to 
integrate down to the continuum, which will require higher S/N than to obtain the redshift of a quasar QSO. At this 
point we also know less about possible relevant variations in galaxy continua. 

In the short term, a pilot study aimed at simply detecting the BAO feature in the Lya forest would be very 
desirable. The bare minimum requirements are roughly > 30 square degrees with a magnitude limit g > 21. 
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